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The ultimate accuracy as regards controlling a qubit with a coherent field is studied in terms of 
degradation of the fidelity by employing a fully quantum mechanical treatment. While the fidelity 
error accompanied by n/2 pulse control is shown to be inversely proportional to the average photon 
number in a way similar to that revealed by the Gea-Banacloche's results [8] our results show that 
the error depends strongly on the initial state of the qubit. When the initial state of the qubit is in 
the ground state, the error is about 20 times smaller than that of the control started from the exited 
state, no matter how large N is. This dependency is explained in the context of an exact quantum 
mechanical description of the pulse area theorem. By using the result, the error accumulation 
tendency of successive pulse controls is found to be both non-linear and initial state-dependent. 
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I. INTRODUCTION 

Nowadays several two-state systems have been pro- 
posed as candidates of physical qubits, and the means 
of single-qubit control and inter-qubit control are also 
investigated. Amongst them, controlling matter qubits 
with coherent electromagnetic field is of particular in- 
terest such as quantum-dot qubit control pQ, spin qubit 
control 2 , single-atom qubit control [3] , ionic qubit con- 
trol @], and Josephson junction qubit control [5] using 
optical or microwave pulses. In practice, we can use the 
classical field theory which enables us to control a qubit 
perfectly as described by the Bloch rotation according to 
the pulse area of the field [51 [7] . 

However, classical fields exist only as an infinitely 
strong intensity limit. In reality, strong but finite 
strength coherent fields are available that will cause con- 
trol errors because of their lack of infiniteness. 

Gea-Banacloche reported that the fidelity error of 
qubit control is limited by a value inversely proportional 
to the average photon number of the field N [5]. A fi- 
delity error of ^Jy 4 is obtained for n/2 pulse control of 
a qubit and directly reflects the coherent noise inherent 
to the control field |S]. Meanwhile, Ozawa considered a 
quantum mechanical bound of the control field both for 
single qubit rotation and two-qubit control gates, and de- 
rived a universal error bound for the single qubit rotation 
to be 16 (n +1 ) from the uncertainty principle [5] , which is 
a rigorous bound but not tight enough for specific models 
including this case. As regards the coherent qubit con- 
trol, Ozawa's bound is ascribed to the phase error part 



of the coherent noise [10] but cannot give the value of the 
total error. 

In this paper, we formulate the full quantum mechani- 
cal interaction between a pure coherent field and a qubit 
in the general initial state including mixed states. Then 
the error rates will be shown to depend on both the field 
and the initial state of the qubit [11]. We break down 
the result into the first correction to the classical pulse 
area theorem by representing the error rates by the or- 
der of 1/N. We also show the successive control of a ir/2 
pulse with a single n pulse and investigate how the er- 
ror accumulation differs from that in the classical case. 
Furthermore, the pulse area theorem is deduced as the 
classical limit of the map by taking the limit N — > oo. 
Here we assume a control field in a good cavity so that 
the Jaynes-Cummings interaction [12] . dominates over 
the interaction of the field mode with the external field 
modes. Our theory is applicable as far as the Rabi fre- 
quency is much larger than the cavity decay constant, 
although we neglect the cavity decay and deal with the 
tendency of the control error for large N cases in this 
paper. In Sec. II, we derive a rigorous solution for the 
qubit dynamics controlled by a coherent field, and ob- 
tain a completely positive trace preserving (CPTP) map 
expression in terms of Kraus operators. In Sec. Ill, we 
derive the general form of the error rate using the CPTP 
map derived in Sec. II. The result will be rendered into 
a simple approximation by the order of 1/N for special 
cases corresponding to the 7r/2 and it pulse controls for 
comparison with previous results. In Sec. IV, successive 
control and error accumulation are examined. Sec. V con- 
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eludes the paper with discussions deriving the maximum 
error rate during the general control of the qubit. 



II. EXACT SOLUTION OF QUBIT TIME 
EVOLUTION INDUCED BY QUANTUM FIELDS 

A. Qubit evolution caused by coherent field 

A control field and a qubit interact during the control 
process and then physically separate. This must intro- 
duce decoherence into each divided state. The coherence 
of the qubit is preserved only when the control field is 
classical, or the infinitely strong coherent field. 

A qubit is assumed to be a two-level system whose 
upper and lower states are assigned as "1" and the lower 
as "0" , respectively. Corresponding state vectors, which 
are denoted respectively as |1) and |0), are written using 
the parameters <f> an d #o (0 < </>, &o < 2ir) as 



\qubit) = cos0 Q |O) + e^sin0 o |l) 



v ' 1 COSfc'o 



where 

and the representation basis is 



e #/2 o 
e" 1 ^ 2 



ID 



|0) 



(1) 
(2) 

(3) 
(4) 



throughout this paper, cos 2 0q and sin 2 6q are the prob- 
abilities of taking the values "0" and "1" , respectively, 
when the qubit is measured in basis Q. 

In our model, the qubit is controlled by a single-mode 
initially-pure state field via a fully quantum mechanical 
Jaynes-Cummings interaction without detuning |12j . Let 
the initial total density operator be 



/ 5(0)=p f (0)®p q (0), 



(5) 



where pi (0), /5 q (0) are density operators for the field and 
the qubit, respectively. The density operator after inter- 
action time t is expressed as 



p(t) = Up(0)U, 



(6) 



where U is the unitary operator of the evolution. 

We assume a non-detuning case, namely, the total 
Hamiltonian is expressed as 



where 



n = hio (a f a + a z /2) + Hi ■ 



1 

-1 



(7) 



is the Pauli matrix representing population inversion and 
a) , a are the creation and annihilation operators of the 
field, T-Lj is the Jaynes-Cummings interaction Hamilto- 
nian, which is naturally deduced from a linear interaction 
model with a bosonic field after taking the rotating wave 
approximation [12 , 



Hi = hg (aa^ 



(8) 



where g is the parameter for the interaction strength de- 
termined by the qubit dipole moment and the local field 
strength, which may be modulated by a cavity, and a + 
and cx_ are the elevation operators for the qubit, which 
are represented by the basis defined by Eq. Q as 



<7_ 



|1)(0| 
|0)(1| 



1 




1 



(9) 



Using the commutation relations for field and spin oper- 
ators, the evolution operator U under interaction picture 
is expressed as a function of k = gt as 

*t*_)] 



U(k) = exp [— in (do"^ 

oo 1 

= E±(J)'(a* + + **-)' 



i=o 



^ (2?)! 

z=o v ; 



2l + 1 i 
+(aa + +a t <r_) 2 '} 



)(a<7+ + a^er_) 



21 + 1 



(2l)\ { i ) 
i=o v ; 

xl 



1 



21 + 1 i 



){(n + l) l aa + + h l a)a-} 



+{a-a + h l + &+a-(n + 1)'}] 



COS Ky/fl + 1 

-ZK(sinc ^^^/n)a) 



-i«;(sinc K\Jh + l)a 
cos K\fh 



(10) 



where h = a^a, 



f smx/ x 
sine x = < ' „ 
1 x = . 



(11) 



First, we assume that the field is in a general pure 
state (not necessarily in a coherent state). Then, p/(0) 
is expressed by photon number expansion as 



Mo) 



oo oo 



rn=0 1=0 



\m){l\Ct 



(12) 



where |m), (Z| are the number state vectors of photon 
numbers m and I, respectively, and C m ,C* G C are the 
coefficients that correspond to the photon numbers m 
and I, respectively. Using Eq. (12 1, the total density 



operator after evolution is written as 

m = u(K)pt(o)®p q {o)uHK) 

oo oo 

= EE^(«)NCmPq(o)qr<i|trt(«). (is) 



m=0 1=0 
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Thus, the density operator of a qubit is obtained in a where the Kraus operator M n is described as 
completely positive trace preservation (CPTP) map as 

oo 

Tr f {p(i)} = 5>|p(*)|n) 
n=0 

oo 

= ^M„p,(0)Mt (14) 

n=0 



oo 

M n = £)<n|#(/s)|m)C n 

771=0 
OO 



E / c>7i,m cos kV m + 1 — i«V?^n+i,mSiiic K\fm \ 

m \ —in-Jm + 15„ m _i_isinc K\Jm + 1 5 nm cosK-v/ro / 

771=0 V ' " 7 

C„ cos K\/n + 1 — iC n +i sin fcyn + 1 
— «C„_i sin K^/n C n cosKV"- 



(15) 



where C_i = is added for the definition of the field 
coefficient. 



Eqs. ( 14 ) and ( 15 ) allow us to calculate any qubit state 



including a fully mixed state evolved by any initially pure 
quantum field. As an example, qubit states after the 7r/2- 
pulse control with a variety of control field types, i.e., a 
number state, a classical state, or a coherent state, are 
shown in Table [J 



B. Pulse Area Theorem as Classical Limit 

When we assume that the control field is a coherent 
state with strong intensity, i.e., the average photon num- 
ber N is very large, the pulse area theorem states that 
the CPTP map ( 14 ) should be well approximated by the 



classical Bloch rotation 



where 



R(<d) 



cos r? sin $ 
— sin i? cos -d 



(16) 



and d = nyN is half of the pulse area of the control 
field. In this section, we examine the N —¥ oo limit 



of (14) to prove the theorem, which is further refined 



in Sec. Ill to obtain an asymptotic expression in terms 
of 1/N. Noting that the field coefficients of coherent 
states satisfy the relation C n +i/C n — yjN/(n + l)e lv 



and |C„| 



c- N N n /n\ for n = 0,1,2, the contri- 



bution of the n photon number component of the Kraus 
operators are written as, 



M n = C n P\<p- 



N 



N 



(17) 



where 



m(x) = 



( 



cos ( §\ 1 + x + i 



^ -VI + a; sin (&y/l + x) cos + x^j j 



(18) 



Choosing <p — tt/2 to eliminate the phase factor in 
Eq.p7|, we obtain from (fl4|), (fl7| and (Il8| that 



lim Trf{/5(t)} = lim |C„| 2 m(a;)/5„(0)77i(a;) 



n=0 



(19) 

In order to calculate (|19|), let us define Hk which corre- 



sponds to the fc-th central moment of n/N as 



2 x k 



71=0 



E< 

71=0 



-N - 



N n fn-N 



n! V N 



Obviously, 



and for k > 2, 



Mo = 1, 
Mi = 0, 



(20) 



(21) 
(22) 



Mfe = e-"^— | nin-Nf- 1 ~N{n-N) k 



71=0 



til 



k-1 



-N 



J\fii — k+1 



k-1 



E 



E 

71=0 



(fc-1)! 
^ i!(fe — 1 — i)! N l 



{(ra + 1- iV)*- 1 - (n- TV)*" 1 } 



(23) 
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Field type 


Pq\t) 


Classical state 


Kii) 

Ideal Control 


Number state 


\ 

Qubit become 


( i o\ 

s completely mixed. 


Coherent state 


V V^ cos ^\^) sin (f\ 


/n + l\ 

/ N ; 


V ^cos(f/f)sin(f V ^I) \ 
cos 2 (fyf) / 



TABLE I: Qubit states controlled from the ground state |0) to an evenly occupied state (|0) + |l})/-\/2 are shown. The number 
state field makes the qubit completely mixed while the classical state (strong coherent state limit) maintains its purity by 
employing a, n/2 pulse. The coherent state with finite strength is expressed by the sum of each photon number component 
where k, is set as ty/{A~^N) to correspond to the pulse area tv/2 . The coherent sate ty/2 pulse cannot give the ideal control as 
far as TV is finite. 



Using ([2lJ)-<[23J), we find that 

1 



Ufc < — , for k > 1. 
N 



(24) 



Eqs. (|2_Lj), (|22j) and (|24|) mean that only the 0-th or- 
der term of x contributes to r.h.s. of ( 19 ) when 
rh(x)p q (0)rh(xY is expanded as a power series of x. 
Therefore, the sum in d 1 9h is carried out as 



JV 



lim Ti-{{p(t)} = lim m(0)p 9 (0)m(0) 



(25) 



which is turned out to be a quantum mechanical proof of 
the pulse area theorem. 



imperfect. To estimate the imperfection of the control, 
we define the error rate by using the fidelity T as below 
after the preceding study [5]. 



P 



l-T 2 



(26) 



III. FIDELITY OF tt/2-PULSE CONTROL WITH 
FIELD WITH FINITE STRENGTH 

When the quantum nature of the control field is not 
negligible, the control of rotation from 9q to 9q + 9 will be 



where \^) is the state orthogonal to the target state|?/>). 
Using Eqs. ( 14 1, ( 15 ) and ( fl8] ) , the error rate for the gen- 
eral control from R{8 )\0) to R(9 +9)\0) can be described 

as 



oo 

P(9 ,9) = (^|^M n p 9 (0)Mt|^) 
n=l 

oo 

= (l\R(-9 -9)Y / MnR(9o)\0)(0\R(-9 )MlR(9 + 9)\l) 

n=0 

GO 

= ^2 I 6 ™' 2 ' sin0 o{ ra11 cos ( o + ™ 01 shi(6> + 9)} 

n=0 

+ cos 9 {m 1Q cos(6» + 6) - m 00 sin(0 o +6)}^ (27) 
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where m 1J = (i\m(x)\j) as shown in (18). Focusing on the ir/2 pulse controls starting from |0) and |1), the error rates 
are obtained as below. 



P + = P(7r/2,7r/4) 



JV 



sin (Ky/n) + cos 2 (k-\/™ + 1) 
2 



1 °° In 
- - |C„| 2 W — sm(ny/n) cos(kV« + 1) 



sin(Ky / w) cos(n\/n + 1)| 



P- = P(0,7r/4) 



oo 



sin 2 («Vn + 1) + cos 2 (/c^/n) / JV 



n=0 



— sm(n\/n + 1) cos(/t\/??)| 



71=0 



^ ^ / JV 

- - V" |C„| 2 \/ — — am(K\/n + 1) cos(ka/™) 



1 00 r^7 

2 - E l C '»| 2 \/7^ sin ( K V") cos(kV« - 1) 

71=1 " 



(28) 



(29) 



For a fixed $(= «V JV), Eqs. (28) and (29) can be written 
as a single formula as 



oo 



n-N 
JV 



(30) 



where 



= Vl + x sin(rJVl + s) cos f 0^ 1 + x ± ^ j 

(31) 

In order to obtain P ± , let us expand f ± (x) in powers of 
x(= n N ) as 



From (|20j, ([21] and 

= iNT 1 , 



(32) 

(33) 
(34) 



derived from ( 23 ) , we have the expression 



n=0 



n - JV 
JV 



Since, the expansion of (31) gives 



/o 



(36) 



ft = T^sin 2 2. + (JV- 2 ), 

8 16 



we obtain from ( 30 ) and ( 35 1 as 
1 



P 



± 



1 - sin 2-d 

~2 + 16JV 
±40(1 - cos 20)} + 0(N~ 2 ) 



{(l + 4ir)sin20-2i?cos2i? 



(38) 



Taking large JV limit of ( 38 ) as 

. 1 - sin 2r? 
lim P ± = 

JV->oo 2 



(39) 



■d = 7r/4(= 0) is found to be the best for 7r/2 pulse control 
at the classical limit. 

For finite JV, ■& — 7r/4 gives 



(^±2) 2 
64JV 



0(JV" 



(40) 



However, (40) may not give the minimum error rates, 



because 7r/4is no longer the optimal value for finite JV. 
The optimal $ to make P ± minimum are obtained by 
solving 



dP ± 



= 0, 



(41) 



as 



1 /3?r 



JV V 32 



± 



16 



')+0(JV- 2 ), (42) 



where i!^ is the optimal ?9 for P ± . The differences of 
^ from 7r/4 in (42) contribute to P ± in the JV~ 2 order. 
Therefore, (40) is approved as the expression for mini- 
mum P ± . In other words, 0(JV _1 ) tolerance is allowed 
for setting $ when JV is large enough. 



G 



The above argument is solely aimed at optimizing the 
fidelity, and the choice of 1? in ( [42] ) gives the optimal 
fidelity up to the order of O^/N 7 ). As we will soon 
see, this choice results in a bias in (<t z ) of order 0(1/N). 
Hence, there may be cases where another choice of $ is 
preferred, which suppresses the bias in (cf 2 ) to the order 
of 0(1/N 2 ) while achieving the optimal fidelity up to 
the order of 0{1/N). Instead of P ± , we will consider the 
deviation of the diagonal elements of the qubit density 
operators from the target value 1/2. Using ( 14) and ( 15 1, 
the deviations, A ± , are obtained as 



OO 

y 

n=0 



C n \ 



2$\ 



i±i 

2 



N 



(43) 



When ■& = 7r/4, the leading terms of A 1 * 1 are ( - 2 „^ 7r , 



which vanish by setting = *& , where 



327V 



(1±2). 



(44) 



Although ( 40 ) is an asymptote that the actual error rate 
of the 7r/2 pulse control converges only at N — > oo, Fig. [I] 
shows that it can be a good approximation for even in 
the range N < 100. Fig. [T] also makes it clear how our 
results are compared to the previous results [H [TUj . All 
the asymptotic curves in Fig. [I] are of the order 1/N but 
our curve is not unique but varies greatly up to about 
20 times depending on the initial state of the qubit. The 
curves corresponding to P ± are simply selected as the 
typical condition but their average value happens to be 
the same as the semiclassical result [8]. It is interesting 
that the error rate starting from the ground state P~ is 
so small that the curve stays lower than the curve 



16(JV+1) 

derived from the uncertainty principle. This is not a con- 
tradiction because the curve ,„,}, , ■■, is the lowest limit 

16(iv + l) 

of the worst case for the control [S] and the curve for P + 
in Fig. [I] locates above it apparently. 



IV. ERROR RATE AND DECOHERENCE 
AFTER SUCCESSIVE FLIPPING CONTROLS 

In Sec. Ill, we quantified the error rates in the qubit 
controlled by a coherent field, which depend not only 
on the strength of the control field but also on the the 



initial condition of the qubit. Since (14 1 allows any qubit 



state including a mixed state as its initial state, we can 
calculate the error rates of successive controls. Here, we 
consider a simple case, i.e., successive ir/2 pulse controls 
to the qubit initially in the ground state or in the excited 
state. 

Let PJ denotes the error rate for the successive 7r/2 
pulse controls that intend to make 9 = — s- 7r/4 (it/ 2 
pulse) and 6 — ir/4 — > ir/2 (ir/2 pulse). If we set # = fl^ 
instead of ir/4 for each control, (d z ) vanishes up to 1/N 
order. Therefore, the density operator after the first con- 
trol = —> 7r/4 can be approximated as the statistical 



Average Photon Number: N 




Ozawa's bound" 
7 Initially in down state 



FIG. 1: The solid lines are exact error rates (1 — fidelity 2 ) 
plotted to the average photon number of the field iV on a log 10 
scale for the initially "0"(down) qubit case and the initially 
"1" (up) state qubit case (from bottom to top). Both show 

good asymptotic agreement even for a small N of ~ 100, n 64 ^ 



(Gea-Banacloche's result) 
to bottom. 



' 16(JV+1) 



(Ozawa's result), from top 



mixture up to 0(1 /N) as 

Ti f {p(t)} = {1 - P(0,7r/4)}|7r/4)(7r/4| 
+P(0, 7r/4)|37r/4) (3tt/4| + 0(1/ TV 2 ), (45) 

where |7r/4) = — - is the target state and |37r/4) = 



is its orthogonal state. The total error rate after 



10) - |1) 
V2 

the second 7r/2 pulse control is then written simply as 

P- = P(7r/4 l7 r/4){l-P(0,7r/4)} 

+{1 - P(37r/4,7r/4)}P(0,7r/4). (46) 

We will calculate P d ~ to the order of 0(1/ N). P(0, tt/4) 
was obtained in Eq. ( 40 1 in the previous section as 



P((W4) = P" = + 0(1/^V 2 ). (47) 



P(7r/4,7r/4) is also obtained using ( |1S| and (27) as 

1 oo 

P(7r/4,7r/4) = 9 ^|C„| 



|2 |m 00 + m 01 | 2 



n=0 



oo , 

-y/l + X sin (|Vl + 



647V 



0(l/iV 2 



(48) 



Finally, P d is expressed by the simple sum of the error 
rates of the successive 7r/2 pulse controls up to the order 



7 



of N- 1 as 



Putting the elements mjf in ( 18 ) , Eq. ( 53 1 for the specific 
values of 6q 



0, 7r/2, ±7r/4 becomes 



P 



d - P(0,7r/4) + P(7r/4,7r/4) 



327V . 



(49) 



Similarly, let denotes the error rate for the successive 
tt/2 pulse controls that intend to make 9 = tt/2 — > "Sir /A 
(it/2 pulse) and 6 — 'in /A ->• tt (tt/2 pulse). We find that 



P 



- P(7r/2,7r/4) + P(3tt/4, tt/4) 
= P+ + P- 



32JV 



(50) 



( 49 1 and ( 50 ) show that the error rates in the successive 
tt]2 pulse controls coincide with the semiclassical result 

x 2. This may be understood as that the total 
controls in both cases (|0) O |1)) can be regarded as 
classical flipping and thus the back action from the qubit 
becomes invisible. With the exception of these special 
cases, successive control errors cannot be explained by 
semiclassical analysis. For example, the error rates of 
the successive tt/2 pulse controls starting from \tt/A) and 
|37r/4), denoted respectively as P d and P d 9 , have differ- 
ent values as 



P(0,tt/2) 
P(7r/2,7r/2) 
P(±tt/4, tt/2) 



^|C„| 2 cos 2 (^) ; (54) 

oo 

|C„| 2 cos 2 (kV^TT), (55) 

n=0 
_^ oo 

- ^ |C«| 2 { cos(ny/n) + cos(ny/n + 1) 

± \/N/{n + 1) sm(n\/n + 1) 
^fy^n/N sin(KV / n)} 2 

1 i 00 r 

2 + 2 X] l^"! l cos(«\Ah- r)cos(/tv^) 

n=0 



/ 

- 1 / sm(fcv ^ + 1) sinffc-v/ n) 

V i + n 



~F\I sin(n\/n){ cos(kV n-j-T) 



i(n\/n — T)} 



(56) 



Using x = (n — N)/N again and setting z9 = KviV = 
n/2 + Q(l /N), w e obtain the expansions for the sum- 
mands in ( 54p6 ) as 



cos 2 (kv / "') = - 1 1 + cos (ny/ 1 + x^j | 



16 



(57) 



P/| ~ P(tt/4, tt/4) +P(tt/2, tt/4) ^^J^ (51) 
P d 3 * ~ P(37r/4,7r/4) + P(0,^/4)~ (52) 



! (kV^TT) 



i|l + cos^l + , + lj| 

2 2 

— x 2 + — (x-3a; 2 ) (58) 
16 327V v ; ' v ; 



We can also calculate the error rate for the it pulse from 



(27) as follows. By setting 8 — tt/2, (27) is rewritten for 



tt pulse control for the initial state as 



oo 

P{8 ,tt/2) = i^|C n | 2 {m 00 + m 11 + (rn 01 + m 10 )sin2^ 



cos( K\Jn + 1) cos{K\fn) — A sin( nsfri + 1) sininy/n) 

V 1 + n 

cos(|Jl + 2; +l) cos(^—x 



n=0 

+(m 00 -m n )cos26lo} 2 



(53) 



1 + X . / 7T , , 

- ,, / 3- sm — < / 1 + x , 

1 tt 2 2 7r 2 (x-2a; 2 ) 
~ 1+ 2N + Y X + 87V 



)sin ( — y/l + x 
\2 



(59) 
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FIG. 2: Error rates of the tt and tt/2 controls(6< = tt/2, tt/4, 
respectively) from the initial values 9o = and tt/2 are plotted 
for the target values 8o + 9. Open squares on the two kinked 
lines represent the error rates for the successive tt/2 pulse con- 
trols with the average photon number TV in each control pulse, 
while the open circles are for the single tt pulse controls with 
average photon numbers TV(uppers) and 27V (lowers), all ob- 
tained from the quantum analysis. The solid diamonds repre- 
sent the corresponding error rates obtained from the classical 
analysis. The error rates obtained from quantum and classical 
analysis show good agreement except the kinked points. 



Sin 



TV 

Vl + i sin 



(kv^i){ cos(k-\Ai + 1) — cos(Ky/n — 1)} 



cos + 




2N 327V. 



(60) 



FIG. 3: Error rates of the tt and tt/2 controls(# = tt/2, tt/4, 
respectively) from the initial values do = tt/4 and 3tt/4 are 
plotted for the target values 9o + 9. The legend is the same as 
used in FIG. 2. In this case, successive tt/2 pulses show linear 
error increase even with the quantum analysis. The quantum 
results are totally different from the classical counterparts and 
show large discrepancy between the two initial conditions. 
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tt/4 


tt/2 


3tt/4 


tt/4 


(*- 2) 2 
647V 


(n+2) 2 
647V 


{tt+2) 2 
647V 


(tt-2) 2 
647V 


tt/2 


n 2 
167V 


(tt+2) 2 
167V 


n 2 
167V 


(n-2) 2 
167V 



TABLE II: Error rates P(9q, 9) of 7r/2 and tt pulses for various 
initial states up to the order of TV -1 . Effect of the qubit back- 
action appears in all the cases except those of tt pulse with 
6*o = 0, tt/2, which corresponds to the classical bit flip. 



ignored the terms of 
2 ) in the approximations. 



where we have 
0(x 3 ),0(xiV'~ 1 ),0(iV' ~) m 
Eqs. (|54| )-(56) are carried out in the same manner as 
Eq. ( |35| ) with regarding the expansions |57|-(60| as 
that of f ± in (32). Thus, the error rates in the tt pulse 



controls up to the 1/tV order are obtained as 



P(0,7r/2) 
P(7r/2,7r/2) 
P(±7r/4,7r/2) 



7T 

167V, 

TT 2 

167V, 
167V . 



(61) 
(62) 
(63) 



The resulting error rates are listed in the Table [TT] and 
shown in the Figs. [2] and [3j which show that the values 
of the error rates and error accumulation tendency differ 
greatly depending on the initial 8q of the qubit. 



V. CONCLUSIONS AND DISCUSSION 

We have analyzed the dynamics of a qubit state con- 
trolled by coherent electromagnetic fields with a fully 
quantum Jaynes-Cummings interaction. Although the 
control is assumed to be free from technical imperfec- 
tions such as frequency detuning and interaction time 
fluctuation, the control error inherent in the finitcness of 
the control fields is found to limit the precision of the 
control. The error rate in tt/2 pulse control with the co- 
herent state of the average photon number TV has turned 
out to be v 64A / by the measure defined by ( 26 ) when the 
qubit is initially in the ground state. It is about 9% of 
the previously reported error rate, obtained from a 
semi-classical model 8J . Since this control can be used to 
prepare the initial state of the qubit register, a one order 
of magnitude increase in the accuracy (or 90 % energy 
saving) is good news. When the initial state of the qubit 



9 



is the excited state, the error rate in the same control be- 
comes . The semi-classical model underestimates 

o4iV 

the error rate by 1.9 times in this case. 

We also have shown the error accumulation in controls 
using two successive 7r/2 pulses. In both cases of a qubit 
starting from the ground state and from the excited state, 
the overall error rate coincides with the semi-classical 
result as 2 x \^ except that the accumulation tendency 
is not linear, as shown in Fig. [2j It is also shown that the 
total error rate is smaller than the single n pulse control 
with a coherent field containing N photons on average, 
but is larger than the 7r pulse control with a field with 
2N photons, which uses the same total energy. 

As shown in Fig. |3j the error accumulation of the suc- 
cessive 7r/2 pulses are linear in the cases where the ini- 
tial qubit state vector lies in the XY plane in the Bloch 
sphere. The successive 7r/2 pulses yield the same error 
rate with the single 7r pulse if the total energies in the 
control fields are the same. But the error rates can differ 
by as much as 20 times, depending on the initial state 
of the qubits. This suggests that the error rate takes 
its minimum and maximum values when the midpoint of 
the rotation of the Bloch vector is located at the bottom 
and the top of the sphere, respectively. This conjecture 
is readily confirmed by calculating the general error rate 
P(9o,0) in (27) up to the order 1/N, which can be per- 
formed in the same manner as was done for P(0, 7r/4) 
and P(n/2, n/4) in Sec. Ill and IV. By using the altitude 
angle $ of the midpoint of the rotation instead of 8o, the 
result is expressed as 



Eq. (64) gives <& values that provides the maxi- 



mum(minimum) value of the error rate at 7r/2(0) for the 
domain < 6 < it. The maximum and minimum error 
rates in the 6/2 pulse control are written as 



PmaxW 



and 



AN 



- smf 
AN 



(65) 



(66) 



The value of -P m ax(#) is important because it is used to 
estimate the maximum error during quantum computa- 
tions where qubit states are unknown before control. 

To conclude the paper, we discuss the validity of the 
pulse area theorem. We have proven the theorem rigor- 
ously up to the order 1/N. If the deviation of ■& from 
7r/4 is in the order of iV _1 , the error becomes quantum- 
mechanical limited and solely stems from the finiteness 
of the control field. If we consider the error rate of the 
order of N~ 2 , the minimum value will be attained at 



cos 2$ sin 9) 2 
AN . 



(64) 
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